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The Proportionality between Charge Acceleration and Radiation 
from a Generic Wire Object 


Edmund Miller™ 


Abstract—The Lienard-Wichert potentials show that radiation is caused by charge acceleration. The 
question arises about where charge acceleration occurs on the most basic of antennas, a center-fed, 
perfectly conducting dipole, for which there are two obvious causes. One is the feedpoint exciting 
voltage that sets into motion an outward-propagating charge and current wave at light speed c in the 
medium. A second is at the dipole ends where the outgoing wave is totally reflected producing a change 
in charge speed of 2c. In addition there is the decreasing amplitude of the propagating wave with distance 
due to its partial reflection along the wire. That reflected charge also undergoes a speed change of 2c. 
This is the reason why the decay of current flowing along a straight wire antenna has been attributed 
as being due to radiation. Radiation caused by these and other kinds of charge acceleration due to 
resistive loads, right-angle bends, and radius steps are investigated. These phenomena are examined 
primarily in the time-domain where they are more observably separable in time and space than in the 
frequency domain. The current and charge induced on an impulsively excited wire antenna and its 
broadside radiated E-field are computed using a time-domain, integral-equation model. The computed 
results are used to derive a numerical relationship between the amount of accelerated charge and its 
radiated field. This relationship is denoted as an Acceleration Factor (AF') that is obtained for various 
charge-accelerating features of a generic wire object are normalized to that of the exciting source for 
comparison with their respect speed changes. 


1. INTRODUCTION 


The Lienard-Wichert potentials show that charge acceleration causes radiation [1]. Thus, in order to 
develop a quantitative understanding of radiation from a perfect-electric conductor (PEC) excited as 
an antenna, it’s necessary to determine where and how much charge acceleration occurs on it. It is well 
accepted that a dipole antenna radiates from its feedpoint and ends, but there is some disagreement 
about whether it radiates elsewhere [2,3]. Numerous authors [4-17] associate the current decay down 
a long-wire antenna as being due to radiation. This seems plausible since the current decay means 
that the corresponding longitudinal Poynting’s vector then carries less power or energy with increasing 
distance [17]. But if current decay causes radiation this means that charge acceleration is somehow 
involved. Charge acceleration occurs because some charge reflects in a direction opposite to that of the 
propagating wave, similar to what happens at the end of the wire, an observation made by only one 
of the references above [16]. It is clear that wherever charge reflection occurs on an antenna radiation 
results. Feedpoint and various reflection-caused radiation mechanisms are explored quantitatively in 
the following. 
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2. CHARGE REFLECTION AND RADIATION 


Detailed attention in this article is primarily focused on a PEC straight wire for simplicity. However, 
it is illuminating to consider conceptually a more general wire object as sketched in Figure 1. The 
radiation due to charge acceleration associated with each highlighted feature denoted by alphabetical 
letter is briefly listed below. 


(i) A spatially varying wave impedance that produces charge reflection along a straight section of wire, 
or what is called here “propagation radiation”. 

(ii) Any electric field applied to the object in Figure 1, whether from some external excitation such as 
an incident plane wave or from a localized source such as the feedpoint of an antenna that will cause 
a current to flow in response. Setting charge into motion produces accompanying acceleration that 
causes radiation. 

(iii) Charge reflection at an open end of a wire is another cause of radiation. It’s worth noting that 
reflection produces a charge speed change of 2c whereas a source field as at B causes a change of c. 

(iv) Charge reflection from a lumped impedance. 

(v) A directional change in charge velocity at a sharp bend or smoother curvature produces charge 
reflection and charge acceleration that causes radiation. (Velocity is used here when a directional 
effect occurs, but when only straight-line motion is involved the term used is speed). 


(vi) Charge reflection at a change in wire radius. 
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Figure 1. A generic wire object exhibiting various features that cause radiation. 


That charge acceleration causing radiation can be demonstrated graphically using the E-field kink 
model of radiation [18]. This is shown analytically by the Lienard-Wiechert potentials [1] from which 
the radiated component of the electric field due to a charge gq can be written as 


sts {rx |(e- 2) « By b eae dex f(r) xa)|) ay 
where s = r-—(u-r)/c. (1b) 


Equation (1) yields the electric field E of a charge q moving with velocity u and acceleration a along a 
prescribed path in space. There is no radiation unless a is nonzero. Applications of Equation (1) include 
accelerator design and analysis of electron lasers. Charge acceleration (1) is not explicitly included in 
a computer model such as TWTD (Thin-Wire Time Domain) [19] that was developed for modeling a 
PEC. But the effect of acceleration can be derived from the results TWTD provides as will be shown 
below. It is the primary goal of the following discussion to develop a quantitative relationship between E 
and q for the various kinds radiation that can occur on the type of PEC object in Figure 1. This is done 
by conducting computer experiments using two well-validated computer codes, TWTD for time-domain 
results and NEC (Numerical Electromagnetics Code) [20] for frequency-domain results. 

Radiation due to charge acceleration is caused by a source field, and several kinds of charge reflection 
are evaluated in the following. The goal is to develop from computed TWTD results a quantitative 
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estimate of the relationship between the accelerated charge and EF), the broadside, electric field parallel 
to a straight section of a wire that it radiates. The proportionality between the numerically computed 
charge and radiated electric field values from TWTD is expressed here as an Acceleration Factor (AF). 
Thus, the absolute value of || can be written as 


|B\| = AF x |Qel, (2a) 
so that AF = |E)|/|Qe| (2b) 


with Q the total accelerated charge density on the wire associated with a particular |E\|._ Absolute 
values are used here for Qc and FE) since charge pulses of opposite signs simultaneously accelerated in 
opposite directions can produce additive fields of either sign. For the case of source, propagation and 
end-reflection radiation (A, B and C of Figure 1) the pulse shapes from TWTD are used from (2b) 
to illustrate the temporal correlation between Ey and Qc. Otherwise, the numerical AF’ values are 
obtained using the corresponding peak values of these quantities to show the effect of some parameter 
variation such as length, radius or resistance loads. The absolute AF values obtained from the various 
cases investigated are not themselves of primary interest here. It is rather to establish to what extent Fj 
and Qc are proportional in time or to a particular parameter. Also of interest are their values relative 
to that of an exciting source. This is because the source causes a change of c in the charge speed while 
any reflection results in a 2c speed change, a phenomenon that is illustrated further below. 

Unless stated otherwise, a wire radius of 107? m is used throughout together with a segment length, 
Az, in the frequency domain of 0.01 wavelengths and in the time domain of 0.01 m, with the time step 
At = Az/c. Some other materials related to the discussion that follows can be found in References [21] 
and [22]. The excitation is a 1-V peak Gaussian pulse for all time-domain results presented in the 
following. 


3. CURRENT DECAY AND PROPAGATION RADIATION IN THE FREQUENCY 
DOMAIN 


The current on an infinite PEC wire antenna varies with distance z from its feedpoint as [8] 


- oikz Oni 
T(z) = me log ———_—_ 4 (3a) 
le 2Cw + log (ie + [(kz)? +e-77] ++ i5n) 
where Cy = log(1/ka) — y (3b) 
and y = 0.57721566 (Euler’s constant) (3c) 


with & the wavenumber of free space and a the wire radius. 

The envelope of maximum current magnitude from Eq. (3) is plotted in Figure 2 normalized to its 
value at 0.16 wavelengths from the feedpoint for ka = 0.01. It decays monotonically with increasing 
distance from the feedpoint due to propagation radiation (A in Figure 1). By contrast, the outward- 
flowing current on an infinite biconical antenna does not decay with distance, as given by [23] 


I(r) = Ine~**", (4) 


with J, the input current. Because the wave admittance of a conical structure is independent of r 
there is no reflected current and thus no current decay. Consequently, radiation from the infinite bicone 
originates only from its feedpoint where the tips of the two bicones are joined. 

The on-surface longitudinal Poynting vector of a 10-wavelength, center-fed dipole consequently 
also decays with distance as shown in Figure 3 [17]. This is obtained from the dipole’s radial electric 
field (charge density) and azimuthal magnetic field (current density) using NEC. The power flow is 
somewhat oscillatory while declining monotonically with increasing distance from the dipole feedpoint. 
The significance of the current reflection and consequent decrease in power or energy flow is explored 
next. 
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Figure 2. The envelope of maximum current magnitude from Shen et al. [8]. 
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Figure 3. The on-surface, longitudinal Poynting vector obtained from NEC for a center-fed, 10- 
wavelength dipole [17]. 


4. A PREVIEW OF RADIATION AND CURRENT DECAY FOR A STRAIGHT 
WIRE 


Before developing AF estimates for the various cases to be considered, it is instructive to examine 
some initial results for long straight wires. The current J and charge-density Q times light speed c are 
shown at several time steps in Figure 4 on a dipole 599 segments long excited at its center by a 1-V 
Gaussian pulse. It is clear that the positive J and Qc_pulses_are numerically equal on the right-hand 
side of the dipole (i.e., Qc = I) but of opposite sign on the left-half side. This is because a positive 
charge moving to the right produces a positive current, as also does a negative charge moving in the 
opposite direction to the left. Their numerical equality implies that the current and charge carry the 
same amount of energy, an effect that is demonstrated explicitly in Section 12. It’s also clear that 
their amplitudes decay monotonically with increasing distance from the feedpoint. As for the frequency 
domain, this decay exhibits the effect of partial charge reflection. This reflection produces the smaller 
negative current components trailing the positive current pulses on Figure 4. 

Observe that on the impulsively excited, center-fed dipole used to obtain the results of Figure 4, 
the partial reflection of the rightward-moving positive-charge pulse and the partial reflection of the 
leftward-moving negative-charge pulse cause acceleration of both reflected pulses. This acceleration 
produces kinks in their electric field lines in accordance with the E-field kink model [18], to produce 
1/R outward-propagating electric fields as kinks in their respective field lines. Since the charges are of 
opposite sign and accelerated in opposite directions their radiation fields have the same polarity and 
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Figure 4. Charge density Q times light-speed c and the current I for a 599-segment wire excited at its 
center by a Gaussian voltage pulse at several time steps. 
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Figure 5. Incident (black) and reflected (red) (a) Qc pulses near the end and (b) a distance-time plot 
for Qc pulses on a 599-segment straight wire. 


are thus additive. 

The Qc and current pulses in Figure 4 show somewhat indirectly that they propagate at the speed 
of light. The plot in Figure 5 shows more explicitly that this is the case for excitation of the 599- 
segment wire at segment 49 near its left end. Plotted there at time steps 551 and 625 are the initially 
rightward-propagating Qc pulse and its leftward-propagating reflection from the wire end in (a). The 
total distance traveled by the incident and reflected Qc pulses is plotted as a function of time in (b). 
A best-fit straight line to this plot reveals a propagation speed of 2.9868 x 10° m/sec, a value that is 
99.56% of the specified light speed of 3 x 10° in TWTD. Considering the approximations inherent in 
a discretized model, this small difference seems reasonable. The arrow indicates the time at which the 
Qc pulse reaches the wire’s end. Even though there is a slight slowing in the pulse speed during end 
reflection (see Section 12) it is not enough to be discernable in Figure 5(b). 

The radiated, broadside E-field for the time between turning on of a center, impulsive excitation 
and the first end reflection of the outgoing current and charge pulses for various wire lengths is plotted 
in Figure 6. Observe that the first E-field pulse occurs as the exciting source initially sets current 
and charge into motion. The second pulse is produced by the outward-propagating current and charge 
pulses, such as shown in Figure 4, reflecting from the ends of the dipole. Connecting these two radiation 
pulses is a continuous, decreasing radiation field due to reflection of the outward-propagating charge 
pulses, denoted here as propagation radiation. 

Note that propagation radiation from a partially reflected charge pulse propagating along the wire 
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Figure 6. Magnitudes of the broadside radiated E-fields versus time step for impulsively excited, 
center-fed straight wires of various lengths. 


is qualitatively the same as the end-reflection radiation that occurs when the entire pulse reflects from 
the open end of the wire. In both cases the change in speed of the reflected charge is 2c. However, 
observe that the propagation radiation that results is less than end-reflection radiation because less 
charge is reflected in the propagating phase. Also, while their total speed changes are the same, their 
associated accelerations are not necessarily so. This is possible if the time over which the speed changes 
differ because their respective reflection processes differ, something that is discussed in connection with 
end refection in Section 12. That the source-accelerated charge undergoes a speed change of only c 
while the end-reflection change is 2c accounts for their respective radiated peak fields of the latter in 
Figure 6 being almost twice as much. 


5. COMPARING CURRENT DECAY ON A LONG WIRE AND CIRCULAR LOOP 


Current decay over a greater distance and time duration is shown in Figure 7(a). A sequence of current 
pulses at 100 time-step intervals is plotted for a straight wire 1199 segments in length excited near its 
left end at segment 49 by a 1-V Gaussian voltage pulse. Also plotted in Figure 7(a) for comparison 
are equivalent results for a circular loop 1199 segments in circumference. For clarity both plots include 
only the positive current pulses propagating away from the feedpoint to the right for the dipole and 
clockwise for the loop. The speed of the loop-current peaks is ~ 3.0168 x 10°, or somewhat faster than 
that on the straight wire, as can be seen in Figure 7(a). This is apparently because the straight-line 
distance between two points on an arc is less than that on a straight wire of the same length so that the 
interaction fields have a shorter distance to travel. This effect also accounts for the gradual spreading 
out of the leading edge of the loop current. 

The normalized envelops of the peaks in the current-pulse amplitudes obtained from the results of 
Figure 7(a) are presented in Figure 7(b). Note the similarity of the time-domain dipole current decay 
with the frequency-domain plot of Figure 2 above. Current-pulse decay on the circular loop exceeds 
that on the straight wire because the loop curvature produces increased charge reflection. Subsequent 
attention is limited hereafter to the straight wire because its broadside far field is more easily related to 
its propagating current and charge due to the more straightforward time-delay differences from various 
parts of the straight wire. 


6. SOURCE-REGION RADIATION 


The broadside E-field magnitude caused by the feedpoint Qc pulse is plotted in Figure 8 where an 
acceleration factor AF's of 29.4 is obtained, a result that is independent of wire length as demonstrated 
in Figure 6. This independence of wire length occurs so long as the time duration of the excitation is 
less than the propagation time along the wire. On the other hand the exciting-source radiation pulse 
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Figure 7. The rightward-propagating positive current pulses on a 1,199-segment straight wire (blue) 
and clockwise on a circular loop (red) at (a) intervals of 100 time steps excited by a Gaussian voltage 
pulse, (b) whose envelopes are plotted as a function of distance from segment 100. 


is dependent on wire radius as demonstrated in Figure 9. The magnitude of the peak radiation E-field 
pulse for various wire radii is plotted as a function of the charge-pulse peak produced at the exciting 
source. The average source-acceleration factor for the variable-radius case, AF'gp, is 29.2. This result 
confirms the linear relationship expected from the Lienard-Wiechert potentials that the radiated field 
is proportional to the charge. 


7. PROPAGATION RADIATION 


Reflection of a propagating charge pulse in frequency-domain results is exhibited as a decay of the 
outward propagating current of Figure 2 above. In the time-domain results of Figures 4, 6 and 7 it is 
exhibited directly in the decreasing pulse amplitudes with distance and indirectly as a negative trailing 
component of the outward-propagating positive current pulses in Figure 4. 

The reflection of the current and charge demonstrates that the wave impedance presented by the 
wire does not match that of the propagating pulses, as discussed in [9]. The integrated charge in the 
rightward-propagating Qc pulse is plotted as a function of time in Figure 10(a) for time steps T and 
T +1. The difference between the two plots of Figure 10(a) is shown as the amount of charge reflected 
per time step in Figure 10(b). A comparison of the broadside propagation-radiation E-field with the 
reflected Qc charge of a 599-segment wire multiplied by an AF'p of 33.1 is plotted in Figure 11. 
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Figure 8. Time plots of the magnitudes of the broadside E-field and the source-induced charge of the 
599-segment dipole as a function of time to exhibit their proportionality with the a source acceleration 
factor AF'gs of 29.4. 
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Figure 9. The magnitude of the peak broadside electric field due to the exciting-source charge as the 
dipole wire radius is varied over 7 orders of magnitude to exhibit their proportionality with an average 
acceleration factor AF'gp of 29.2. 
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Figure 10. (a) The positive space integrals of Qc at successive time steps on the right-hand half of 
the 599-segment dipole and (b) the difference between successive positive Qc pulses on the right-hand 
half of the dipole in (a) as a function of time. 
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Figure 11. Time plots of the magnitude of the broadside E-field and the reflected charge on a 599- 
segment dipole as a function of time to exhibit their proportionality with an acceleration factor AF'p 
of 33.1. 


8. REFLECTION RADIATION FROM THE END OF A WIRE 


The reflection radiation from a wire end is plotted for a 599-segment wire as a function of time step in 
Figure 12 where the AF'z is 53.5. This is 1.83 times that of the average exciting-source radiation, AF's 
of Figure 9. This ratio reflects the fact that the source-accelerated charge undergoes a speed change 
of c while the end-reflection charge experiences a speed change of 2c. The fact that this ratio is not 
2 apparently demonstrates that end-reflection acceleration is not quite as abrupt as the acceleration 
due to the exciting source. This point is further discussed in Section 12 in connection with the energy 
carried by the current and charge. The end-reflection radiation-field peak was shown to decrease with 
increasing dipole length in Figure 6, due to the decreasing amount of charge that reaches its end due 
to charge reflection along the wire. This is demonstrated quantitatively in Figure 13. The results for 
the propagation, source and end-reflection radiation are illustrated together in the combined plots of 
Figure 14. 
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a 599-segment dipole as a function of time to exhibit their proportionality with an acceleration factor 
AF f of 53.5. 


9. REFLECTION RADIATION FROM RESISTIVE LOADS 


The magnitude of the peak broadside radiated E-field caused by charge reflection from a lumped, 
resistive load is presented in Figure 15 for a wire 199 segments long. A Gaussian voltage pulse is 
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Figure 13. The decrease in the peak broadside electric field caused by a decreasing end-reflected Qc 
pulse with the wire length a parameter. 
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Figure 14. Combined plots of the propagation, source and end-reflected broadside E-field magnitudes 
and the respective time varying accelerated Qc charge pulses that produce them on a (a) linear plot 
and (b) log plot. 


applied at segment 75 with the load located at segment 100. As for the various similar plots above, the 
radiation field is proportional to the reflected charge with an average acceleration factor AF’ for this 
case of 34.1. 
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10. RADIATION FROM DIRECTIONAL REFLECTION AT A 90-DEG BEND 


Reflection radiation caused by a 90-deg bend is plotted with the wire radius a parameter in Figure 16 
for a 199-segment wire excited impulsively at segment 75 with a bend at segment 100. The acceleration 
factor, AF'p, is found to be 30.3. Only a 90-deg bend is examined here as this results in radiation from 
the bent portion of the wire being orthogonal to that from the other wire half. If the bend were different 
from 90-deg, identifying what part of the radiated E-field comes from the two halves of the wire would 
be less clear. The plot of Figure 16 continues the relationship established in the previous results that 
show an essentially linear proportionality between the peaks of the reflected Qc and broadside radiated 
E-field pulses. 
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Figure 15. The peak broadside E-Fields caused by charge reflection from lumped, resistive loads of 
various values. 
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Figure 16. The peak magnitudes of the broadside E-field and the reflected Qc pulse from a right-angle 
bend for an impulsively excited wire 199 segments long with wire radius a parameter. 


11. REFLECTION RADIATION FROM A STEP IN THE WIRE RADIUS 


The remaining example for a radiation feature from the generic object of Figure 1 is a step change in 
wire radius. But in contrast with previous results where the reflected Qc pulses are clearly defined, 
this was not found to be the case for the step-radius wire. Instead, the step effect is more clearly 
demonstrated in a frequency-domain radiation pattern such as shown in Figure 17. The radiation 
pattern of a sinusoidal current filament (SCF) is plotted for comparison together with the patterns of 
two center-fed, 200-segment dipoles modeled using NEC. Each dipole is 10 wavelengths long with a 
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Figure 17. The radiation patterns for a 10-wavelength sinusoidal current filament and two NEC dipoles 
having radius steps at 1/4 and 3/4 of their lengths with each having a 1-A maximum current. 


center section of 100 segments and a radius of 10~° wavelengths. The end sections are 50 segments long 
each and have radii of 10~* and 10~° wavelengths respectively. The current magnitudes for each are 
1-A. 

The SCF and the step-dipole patterns in Figure 17 are similar with respect to their main lobes. 
The effect of the radius steps is to cause 5 smaller lobes to appear at two-lobe intervals between the 
larger lobes beginning at broadside, an effect that increases with an increasing in step radius. These 
results provide indirect evidence for the radiation caused by a step in the wire radius. 


12. AN ENERGY MEASURE THAT DEMONSTRATES RADIATION LOSS 


One way to evaluate the energy radiated by an impulsively excited object is to integrate the far fields 
over an enclosing sphere as a function of time, something that can be computationally expensive. An 
alternative demonstrated here is to instead integrate the boundary sources over the object at each time 
step of the numerical solution [24]. 

This is done here for a wire dipole impulsively excited at its center using 


L 
W,(t)= | I’ («-t)dzx (5a) 
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-£+A z 
7 Q? (x, t)dx + | Q? (x, Ha (5b) 
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where W7(t) and Wg-(t) are the current (or magnetic) and charge (or electric) energy measures at time 
t and A is the segment length in the numerical model. The two terms in the bracketed portion of 
Equation (5b) account for the fact that the propagation speed of the charge pulse going to zero at the 
ends of the wire cannot happen instantaneously. This fact is approximated for by assuming that the 
squared speed falls linearly to zero over the last segment of the numerical model. The c?/3 is an average 


that comes from r J 
1 cx) 2 Cx 
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Observe that if the 1/3 is not used in Equation (5) the total energy no longer decreases monotonically, 
but rises at end reflection, a non-physical result. This result in Eq. (6) also offers some indication about 
how the acceleration term in the Lienard-Wiechert potentials might be influenced by the particular 
details of the radiation features discussed in connection with Figure 1. 
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Figure 18. The current, charge and total energy measures as a function of time for an impulsively 
excited straight wire 199 segments long. 


Results for these two quantities and their sum are presented for a center-fed dipole 199 segments 
long in Figure 18. After the startup transient, W;(t) and Wg-(t) are equal until the J and Qc pulses 
reach the end of the wire where Wg-(t) becomes predominant as the charge accumulates there and the 
current goes to zero. Upon reflection, their equality is reestablished until the oppositely propagating 
pulses meet at the wire center where W;(t) becomes predominant since the current is additive while the 
charge pulses cancel. After each end reflection, the total energy exhibits a marked decline due to the 
end-reflection radiation. This behavior continues periodically until the wire returns to charge neutrality. 
As a result, the total propagating energy decreases monotonically due to the effects of end-reflection 
and propagation radiation. 


13. SUMMARY OF THE DERIVED ACCELERATION FACTORS 


The various acceleration factors developed from the above results are summarized for reference in 
Table 1. 


Table 1. 


ACCELERATION NORMALIZED 
RADIATION TYPE FACTOR AF 


Soares: Variable Wire Rad: AFsx [| 292 ~——SC«dCC“‘i Cd 


14. CONCLUDING COMMENTS 


The results presented here were developed using the computer model TWTD [19] to estimate a 
quantitative relationship between an accelerated charge pulse on a PEC wire and the radiation that 
acceleration produces. The estimates developed for what is denoted here as an acceleration factor (AF) 
cannot be considered numerically precise. This is not because the TWTD results are uncertain, as 
TWTD has been extensively validated. It is instead because there is no clear-cut way for choosing 
at what time a reflected Qc pulse should be “measured” for comparison with the time variation of 
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its radiated E-field. This is true whether the pulse shapes themselves are used or when it is only 
their peak amplitudes are needed to estimate an AF. Thus, the AF’ values given throughout the text 
and summarized in Table 1, being derived from a computational “experiment”, must be considered 
as approximations. Also an issue in this respect is the fact that the actual acceleration depends on 
the specific nature of the physical feature that determines it. This is demonstrated, for example, in 
the values of AF'g and AF’s that might be expected to be in a ratio of 2 : 1 given their respective 
speed-change ratio of 2c/c. The numerical value actually obtained is somewhat less at 1.83. In spite 
of the limitations involved, the numerical results obtained do provide a quantitative insight into the 
relative effects of various radiating features on a generic wire object and do demonstrate a consistent 
proportionality between the radiated E-field and the charge that produces it. 
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